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SUMMARY 


A general integral form of the houndary-layer equation is derived 
from the Prandtl partial -differential houndary-layer equation. The 
generf^ Integral equation, valid for either J.aminar or turbulent incom- 
pressihle houndary-layer flow, contains the Von Karman momentum equation, 
the kinetic-energy equation, and the Loitsianskii equation as special 
cases. 


In an attempt to obtain a practical method for the calculation of 
the development of the .turbulent boundary layer, use is made of the 
e:iqperimental finding that all the velocity profiles of the turbulent 
boundary layer form essentially a single -parameter family. The general 
equation is thereby changed to a simpler one from which an equation for 
the space rate of change of the shape parameter of the turbulent boundary 
layer can be obtained. 

The resulting equation for the space rate of change of the velocity- 
profile parameter is restricted by the assumption that the velocity 
profiles of the turbulent boundary layer can be approximated by power 
profiles. Two of the resulting equations are used to calciilate the dis- 
tribution of the profile shape parameter over an airfoil for one experi- 
mentally determined pressure distribution. Althou^ different assmp- 
tlons were tried for the shearing stress across the boundary layer, the 
calculated distribution of the profile shape parameter did not agree 
exactly with the experimental distribution. 

’ An examination is made of the effect of using the experimentally 
deteimined single-parameter family of velocity profiles Instead of the 
power profiles on certain functions that occur in the equation for the 
space rate of change of the velocliy-profile parameter. One calculation 
of the distribution of the profile shape parameter over an airfoil is 
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also made for the easperlmentally determined pressure distribution by 
using the single -parameter family of velocity profiles found from experi- 
ment. A comparison of the results with those of a calculation made with 
the same assumptions except for the vise of power profiles shows some 
difference near the separation point. It is believed, however, that the 
apparent lack of reliability of the specific equations used to make the 
calculations is caused mainly by the lack of precise knowledge concerning 
the svirface shear and the distribution of liie shearing stress across the 
turbulent boundary layer. The present analysis emphasizes the need for 
Information concerning the shearing stresses in turbulent boundary layers. 


BSTRODUCTION 


An outstanding problem in aerodynamic theory is to calctilate -rfiether 
the flow will separate from the surface of a specific body and, if so, 
where the separation will occur. The concept of the boundary layer and 
the equations that describe the flow in it. Introduced by Prandtl (refer- 
ence 1 ) and first worked out in some detail by Blasius (reference 2), 
reduce the problem to solving the Prandtl boundaiy-layer equation when 
the flow is laminar. Because of the mathematical difficulty of solving 
the equation, approximate methods were developed for the ceJ-culation of 
the properties of the laminar boundary layer (reference 3) • In some of 
these methods, for exanqple, the Pohlhausen method (reference 3) and the 
Wieghardt method (reference 4), a functional form is chosen for the veloc- 
ity distribution through the boundary layer and is combined with either 
the Von Karman momentum equation alone (reference 5) or with both the 
Von Karman momentum equation and the kinetic -energy equation (reference 4) . 
The result is the replacement of the Prandtl partial -differential equation 
by one ordinary differential equation in the Pohlhausen method and by two 
ordinary differential equations in the Wieghardt method. A solution of 
the ordinary differential, equation or equations provides the boundary- 
layer velocity profiles along the body. These and other approximate 
methods that use only the Von Kermm momentum equation, or the momentum 
and kinetic-energy equations together, do not satisfy exactly the Prandtl 
boundary-layer eqv^ation. 

Because the flow in the boundary layer is more often turbulent than 
laminar in cases encountered in engineering, the problem of calculating 
the separation point is of even more importance for turbulent than for 
laminar boundary layers. In spite of the importance of the problem, 
however, less progress has been made in the development of methods for 
the calculation of the behavior of turbiilent boundary layers than for 
laminar boxmdary layers. The lack of progress stems from the absence of 
an explicit independent equation for tl^ shearing stress that is accurate 
enough to lead to a description of the flow when used with the Prandtl 
equation. 
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The main attempts to obtain methods for the calculation of the 
behavior of the incompressible turbulent boundary layer in the presence 
of pressure gradients are those of references 6 to 12. The results of 
these attempts are xinsatlsfactory either because the assumptions upon 
which they rest are incorrect or because the equations used to make cal- 
cTilations were not derived from the boundary-layer equations. 

The analysis of reference 6 is based on the assumption that the 
velocity profile is a single -valued function of the ratio of the pressure 
gradient to the skin friction, an assun^tion shown to be Incorrect by 
later investigators '{for example,' see reference 12).. In the analyses of 
references "J , 10, 11, and 12 the momentum equation is used, together with 
an auxiliary equation, to calculate the distribution of velocity profiles 
over a surface. In each of these four methods the auxiliary equation is 
not derived from the boundary- layer equations but is empirical. 

In reference 8, the equation that gives the variation of the mixing 
length across a pipe (reference 3) was used to calculate the velocity 
profiles. The' fact that the mixing-length distribution across the bound- 
ary layer is not the same as across pipes is shown in references 13 to 15 . 

Reference 9 does not provide a method for the calculation of the 
distribution of turbulent velocity profiles along a surface. It does, 
however, suggest that separation of the turbulent boundary layer always 
occtura when the numerical value of the nondlmensional pressure gradient 
reaches an empirical constant. 

The purpose of the present investigation is to begin with the 
boundary -layer equation for incompressible flow and to proceed as closely 
to a method for the^ calcvilatlon of the behavior of the turbulent boundary 
layer as the present knowledge of the turbulent boundary layer permits. 

At first it might appear that the use of empirical auxiliary equa- 
tions in methods for the calculation of the behavior of turbulent bound- 
ary layers can be avoided by developing a method similar to the Pohlhausen 
method which requires the solution of only the Von Karman momentum equa- 
tion. For turbulent flow, however, in contrast with laminar flow, the 
conditions on the behavior of the veloclly profile at the surface that 
can be obtained from the boundary-layer equation seem to be of little or 
no use for the determination of the shape of the velocity profile across 
the boundary layer. This difference between laminar and turbulent flow 
makes inapplicable the Pohlhausen process in which a type of function is 
chosen to represent the velocity profiles, the function for the velocity 
profiles is combined with the Von Karm^ momentum equation, and the 
resulting ordinary differential equation for the space rate of change of 
the profile shape parameter is solved. 
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■ An auxiliary equation for the calculation of the behavior of the 
turbulent boxindary layer can, however, be obtained from the boundary- 
layer equation by making use of the eacjerimentally verifiable fact 
(references 10, 11, lU, and 15) that all velocity profiles of the 
turb\ilent boundary layer form essentially a single-parameter family of 
curves. In the present analysis the Loitsianskii equation (reference I6) 
is generalized by multiplying the Prandtl boimdary-layer equation not 
only by an arbitrary power of the velocily in the boundary layer but also 
by an arbitrary power of the distance from the surface. The resxilting 
equation is then integrated across the boundary layer and provides a 
general Integral form of the boundary-layer equation, valid for either 
laminar or turbulent flow. This general integral' form of the boundary- 
layer equation reduces to the Loitsianskii equation when the distance 
from the surface is raised to the zeroeth power, to the Von Ka rman momen- 
tum equation when both the distance from the surface ant^ the velocity are 
raised to the zeroeth power, and to the kinetic-energy equation when the 
distance from the surface is raised to the zeroeth power and the velocity 
is redsed to the first power. 

Ifhen use is made of the assimiptlon of a single -parameter family of 
velocity profiles, the general integral form of the boundary-layer equa- 
tion becomes a general equation for the rate of change along the surface 
of the velocity-profile shape parameter. This equation. for the rate of 
change of the velocity-profile shape parameter is the desired auxiliary 
equation. 

The assxnnption of the single-parameter family of velocity profiles 
changes the problem from one of finding a solution of a peirtial- 
differential equation, the Prandtl boundary-layer equation, to one of 
finding a solution of two simultaneous ordinary differential equations, 
the equation for the rate of change of the shape parameter and the 
Von K&rm^ momentum equation. The differential equation for the rate of 
change of the shape parameter, however, ceinnot result in a solution of 
the problem in the present analysis beca'use a knowledge of the shearing 
stress is lacking. In the present analysis 'various assumptions are made 
for the distribution of shearing stress through the boundary layer, and 
the distribution of the shape parameter over the surface of an airfoil 
is then calculated. Because of the arbitrary assumptions for the shear 
distribution and the use of a flat-pla'te skin-friction formula, precise 
agreement between the calculated and experlmen-fcally obtained distributions 
of the shape parameter is not obtained. 

‘The problem of finding the shearing stress in the turbulent boundary 
layer remains . It is believed, ho'wever, that if suitable approximations 
are. found for the shear and sirrface friction, the equations presented 
herein should enable the development of the ■turbulent boundary layer to 
be calculated "with an accuracy sufficient for engineering purposes. 
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The present work was begun while Dr. Lin was temporarily at the 
Langley Laboratory and was continued by correspondence. 

SYMBOLS 


A 

a,b,c 

B 

c 


Fo.Fl 



g = 


jr 

^0 


51 

52 

e 

Ho 


arbitrary positive integer in shear polynomial 


coefficients in polynomial for 



tlfdt 


ejcponent in expression for shear 
reference chord 
functional notation 


derivative of shear polynomial for 1 = 0 

ha 

coefficient of 1 in expression for 


equllibrlvim value of H for o> = 0 
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j = -i- r (1 - f)ay 

gn+1 Jq Jo 


dy 


pi 5-1 N 
/ ^ 
Jo 


J,0;I 


K 


coefficients in polynomial for 


cL^ 


ratio of kinetic-energy thickness to mamentum thickness 

(1/0° ^ 


K- 

dH 


function of H 


(=^ 


M 


= -1— I t(l - f“"^)y^ dy 
gn+1 Jq ' 

= -i- / yii(i _ i“+l)dy 
gn+1 Jq 


m 


eacponent of u in derivation of general equation 




B- =2 

dH 


P 

Pi 


exponent of y in derivation of general equation 
coefficient of cu in equation for 0 

dx 

exponent in equation for power profiles, f = ^P 
static pressure 
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en Jq 



ro 

S 

U 


u 


V 


^0 

X 


y 

5 


6 * 




radius of "body of revolution 

coefficient of 0 in equation for 6 — 

dx 

velocity parallel to surface and at outer edge of 
■boundary layer 

velocity parallel to surface and inside 'boundary layer, 
positive in direction of positive x 

velocity perpendicular to surface and inside boundary 
layer, positive in direction of positive y 

value of V at y = 0 

coordinate parallel to surface, positive in direction 
from leading to trailing edge 

coordinate i>erpendicular to STurface, positive outward 
from surface 

smallest value of y for which the difference U - u is 
negligible 
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Tq dx 



viscosily 

1 = ^ 
° H-l\ 

- I/o' ^ 4 

P 

density 

T 

shearing stress 


surface shearing stress 

pu^' 


II 


U dx 



- , AHALISIS 

Derivation of General Equation 


The general equation is derived for the hody of revolution hecause 
the equation for two-dimensional flow can he obtained from this equation 
hy letting the rad.ius of a transverse section of the hody of revolution 
become infinite. 

The houndary-layer equation of motion for the hody of revolution, 
slIso valid for two-dimensional flow, from reference 3 is 


dx. ^ p dx p 5 y 


( 1 ) 
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After multiplying throu^ by mailng use of the equation of conti 
nuity that is valid in the boundary layer of a body of revolution 
(reference 3 ) 

^+^ + ^-^ = 0 
8 x Sy 1*0 dx 

and noting that 

pD ffi = • ' 

dx dx 


equation ( 1 ) becomes 


+ ^ + 2 i. + u 5u°^^ ^ V ^ 

m+l\Sx 9y Todx/ m+lSx m + 1^ 

u.m/u + i (2) 

\ ax P 

After equation (2) is written in a form in which each term vanishes at 
the outer edge- of the boundary layer, each term of the equation is multi- 
plied by y^ and integrated from y = 0 to y = 6 . The result (see 
appendix A for detailed development) is 


(n + 1)W — + ef— - nl]) + — 2?[N(m + 2) - n(j - M) - L(m + l)] + 
dx 7 U dx 

PB/e aJL 

1 . M): - n« (3) 


Equation ( 3 ) is the general integral form of the boxindary-layer equation. 
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The- Von Karmw mofmentum equation Is obtained froia equation (3) by 
letting m = 0 and n = 0; the equation for kinetic energy Is obtained 
by letting m = 1 and n = 0, and the equation for moment of momenttna 
Is obtained by letting m = 0 and n = 1. 

In the case m = n =. 0, 


Nd 



dy = 0 



or 


L = 



-H 


It can be easily verified that all the integrals, except Q, involved in 
equation ( 3 ) have finite integrands as n approaches 0. The limit nQ, 
however, approaches unity as n apprqaches 0; thus 

= n 


= [yn(i . ® yn L(i . f«-X)ay 


/ y“'^(i - f“*^)ay 


vn 


The first term drops out ,lf n 0 . Then, by taking the limit n -> 0, 
the result Is 


Itm nOe^ = 1 
n— ^0 


Hence, when m = n = 0, equation ( 3 ) becomes 


dx 


e dU/„ 


^ 0 
ro dx 
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Equation ( 4 ) is the momenttmi equation for flow over a body of revolution 
with flow throu^ the surface. For two-d iTtiens ional flow, equation (U) 
becomes 


^+1^(h + 2 ) -!2 = Z ^ 

cbc U dx U 


( 5 ) 


the value 

of 

~ from 

equation 

(4) is substituted into 

equation (3), 

result is 


dx 





- nil) 


^(m(N - 

> 

L) - n[(.: 

r - M) + N(H + 2 )] - L - 

M?- - 


U 

dxl^ 

J 

iq ax '■ 

( J - 

M) + ■ 


N(n + 1}] 




1) 

pb/e 

/ A’ 
'0 


(n + 1)N 


(6) 

PH^i: 

0 

bn 





where 



The assim^tions contained in equation (6) are the usual boundaxy-layer ' 
assun5)tions. Equation (6) is valid, for both laminar and turbulent flow. 


Form of Equation (6) for Single -Parameter Family 
of Velocity Profiles 

Equation (6) is now to be placed in a form valid when the velocity 
profiles form a single-parameter family of curves (f = f(ti,H)). For 
this purpose the term Ii of equation (6) is modified in the following 
manner: 

By definition, 

= J' yii-1^1 _ fm+l) 


/ 

Jo 


^ d(l - f) 


dx 


y 


dy dy 
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Because f depends only on 11 and 


8(1 - f) ^ 8(1 - f) ^ ^ 8(1 - f) ffl 
8x 8ti 8x 8 h dx 


From "the definition of ■q 


8q _ n d0 

8x e dx 


then 



S(1 - f) 
8x 


ay = 


-i“ , S(1 - f) ay + S a(i - f) 

0 dxJo 8tj dx Jo 8H 


or, after an integration hy parts of the first term on the right-hand 
side, the result is 



B(1 - f) 

8x 


ay = .1^ 
0 dx 


y(i - 



f)dy 


ffl r^ay 

^Jo 8H 


The e^qresslon for 


13,9“+^ 


then hecames 


yn(i . f"^l)ay .i| ^ y 3 .f (1 . f«l)ay + 

g^%-l(l.f"*l)(^|ay)ay 
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But 




f yn(i - iJn.+l)f 3y ^ 

'Jo 


W 


n+1 


and 


J y^-l(l - ~J^ (1 - f)d^dy = je^^+l 


Then with 




le 


n+1 


the expression for can he written as 


= i §§(_M + N + j) - 1 as 

0 dx dx 


dH 


C7) 


When the egression for Ii from equation (7) is substituted into equa- 
tion (6) and equation (4) for ~ is used, the following equation is 
obtained: 

® = I - m)] + 


Afn 

pu^L 


n(J - M) - N - (m + 1) ^ di] + 

v>0 ^ _ 


^[n(J - M) - H + n^ 


(8) 


o 
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where — = — _ has heen used. Equation (8) for 0 — is applicable 
'dx an dx dx 

both to tvro-dimensional flow and to flow over a body of revolution. 


In equation (8) all the integrals, except the one involving the 
shear ratio g, are functions of H, m, and n only. For the present 
no restrictive assun^tlons regarding the shear are made. The form of 
the Mnetic-energy equation for a single-parameter family of velocity 
profiles is obtained from eqmtion (8) by placing m = 1 and n = 0 
and dividing by H* (H) = K' (H) . Thus, 


p dH K(H - 1) 0 dU 


0 ^ 

Ia_l 

- 1 \ 

dx K' U dx 

V K. • y 

PU 2 ' 

( K* 1 



The symbol K represents the ratio of the kinetic -energy thickness 
8 

(1 - f2)f ^ -to the momentum thickness. Note that in the derivation 

0 


I 


of equation (9) from equation (8), the assumption of a single-parameter 

VfX 

family of curves is not restricted to the case = 0. 


Restriction of General Equation to Power Profiles 

• The data in figure 1 show that the power profiles defined by f = 
are a good approximation to the "standard" profiles derived by fairing 
experimental data (reference 10) . Equation (8) can be further developed 
by using the assimiptlon that f = ^P. After some fairly lengthy calcula- 
tions (see appendix B), equation (8) becomes 

g _ -kp(p + l)(2p + l)[p(m + 2) + n + 1] £ ^ ^ 
dx pm + n + 1 Udx 


2p[p(m + 2) + n +,,fl.-|( 2 p + 1) + |^(m + 2) + n + l] ^P“+n dO ^ 
pm + n L '•'0 pU^ 

2p(p + l)[p(m + 2 ) + n + l3 ^ 
p(m + IJ + n U 


( 10 ) 
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As a first approximation the asstimption has been made that f = even 

when — 9^ 0. 

U ^ 

The occTirrence of the arhitoary positive integers m nrui n in 
'equations (8) and (lO) requires an explanation. Iii order to determine 
why m and n appear, equation (8) is written in a different form.. By 
making use of the definitions for N, I, J, M, L, and Q and inte- 
grating hy parts where necessary in order' to eliminate terms that con- 
tain the result is 


dH 




-n(H + 1)(J - M) + N(H - m) + L(m '+ l) = 

f dq + f 2 - dq 


p6/0 . 

-n(j - M) - N - (m + 1) / ^ 

Jo ^ 


dq = 


and 


n(j - M) - N + nQ = (m + 1) q“fHi |f ^1 f dq^ 


dq 



Equation (8) then becomeB 



By using the assumption of a alogle-paremeter family of curyes directly in the partial- 
differential equation (l), the following ordinaiy differential equation is obtained: 


dx 




f dTl + 

+ 1 

1 - f2 - (H + 1)|? / 

.n 1 

f dll 

.!2 


f dq - 1 

CM 

[dj\ J( 


U dx 

^ Jc 

) ] 

U 

i^Vc 

) A 


f M . 5£ f g ^ 


( 12 ) 


The concept of a single-parameter family of velocity profiles is consistent with equation (l) 
and with particular functions for tq/pU^, g, and f ■ when the right-hand side of equation (12) 
is independent of t). When the rl^t-hand side of equation (12) is independent of t), the right- 
hand side of equation (ll) is Independent of m and n. Equations (ll) and (12) are then 
Identical . 

JTT 

To obtain an equation for Q — that does not contain either m or n or both, the 

dbc . 

functions Si and f must therefore be such that the right-hand side of equation (12) 

is independent of the solution of the equation for 0 — then provides a solution of equa- 
tion (l). Note that the problem is to find a solution not of equation (l) alone but of 
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equation (l) and the Independent relation for the shearing stress in 
turbulent boundary layers; this relation is at present unknown. 


The nature of the approximation made in the present analysis^ in 

order to obtain a specific equation for — , may be clarified by noting 

dx 

that a specific equation for ^ is obtained from equation (8) by choosing 


the functions Tq/pU^, g, and f and substituting an arbitrary positiye 

integer for m and an arbitrary positive integer for n. The calculated 
distribution of H over a body for arbitrarily chosen functions for 
tq/pU^, g, and f is then consistent with the momentum equation and one 

JTT 

of the Integral equations for — . For exanq>le, if m = 1 and n = 0, 

dx 

both the momentum and the kinetic-energy equations are satisfied but no 
other ones. If m = 0 and n = 1, only the momentum and the moment of 
momentum equations are satisfied. In the present analysis only the 
momentum^ the kinetic-energy^ and the moment of momentum equations - 
equations which have familiar physical meaning - are used. 


As noted previously, equation (U) is independent of m and n if 
the functions Tq/pu^^ and f are such that the right-hand side of 
equation (12) is independent of tj. In this case a solution of equa- 
tion (1) results and the functions Bf aiid S and the calculated 

distribution of H satisfy every particular equation obtainable from 
equation (ll), (lO), or (8) by assigning positive integers to m and n. 

Note that m and n cannot both be made zero in equation (8) 
dN 

because — + nl = 0 for m = n = 0. If m and . n are both zero, 
dH ' 

equation (8) becomes 0 = 0. It is also noted that equations (8) and (?.0) 
are valid both for flow over a body of revolution and for two-disiiensional 
flow. 

For m = 1 and n = 0, equation (lO) leads to the equation for 
kinetic energy 

0 ^ = -h(H - l)(3fl - + (3H - 1 )(h + ^ + 

dx U dx \ ^ Jq bt J pxj2 


(H + 1)(3H - 1) 30 

k U 


(13) 
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where the relation for power profiles 2p + 1 = H has been introd\iced. 
This form of the energy equation can also he ohtained from equation ( 9 ) 
hy noting that, from the definition of K and the equation for power 
profiles, 

• 

K - 2(gP + 1) = 

3p + 1 3H - 1 


A comparison of the values of K ohtained from this formula and Ohtained 
from the standard profiles is given in figure 2. 


The equation of moment of momentum for power profiles is ohtained 
ftom equation (lO) hy letting m = 0 and n = 1; it is 


g ^ ^ H(H + 1 )(h 2 - 1) £ ^ + /g2 1 ) 

dx 2 U dx ’ 


H + (H + 1) / 



3 ^ 


1 °. 

pU2 


(h2 - 1)^ 




In this equation the term involving the shear distrihution may he 
rewritten as follows; 



It then Involves the mean shear inside the boundary layer. 


Attempts to Derive a Relation Governing the Change of 
the Form Parameter 

In most of the recent analyses of the development of a turbulent 
boundary layer, an empirical relation governing the change of the form 
parameter H is usvially introduced. It is clear that equation (lO) 
automatically furnishes such relations if the shear distribution is haown. 
In this section, three attempts are described to establish such a relation. 
These attempts are based on the following simple assunqptlons for the shear 
di stribution : 



NACA TN 2158 


19 


g dpn 

(a) The shear dlstrihution depends only on — , which is equal’ 

■’’O 

to the Pohlhausen parameter multiplied by a factor (reference 3 ) 


(b) The shear is constant across the boundary layer 


(c) The shear distribution depends only on the form parameter of 
the velocity distribution 


The first two assumptions are used either with the energy equation 
in the forms given by equations (9) and (I3) or with equation (l 4 ) for 
the moment of momentum. The last assunqotion is used with equations (13) 
and (1^) jointly. 


(a) Shear distribution depending only on the Pohlhausen parameter . - 
The first assun^jtlon follows the original idea, of the method of Von KarmAn 
and Pohlhausen in using polynomial approximations together with the bound- 
ary conditions obtained by successive differentiation of the equations of 
motion (reference 3)- Fediaevsky (reference 8) appears to have been the 
first to introduce it into the investigation of turbulent boundary layers. 
When the shear stress through the turbulent boundary layer is assumed to 


be a polynomial of fifth degree in S = =• satisfying the following bound- 

o 

ary conditions; * 

« 

at y = 0 


T = tq, 


^ dx * 



at y = 6 



Sy2 


0 


the following expression is obtained; 


g = (1 - + (3 + + 3(2 + x)sd 

The shear distribution g is a function of ^ and where 

■'b U dx 0 Tq 


(15) 



20 


KACA TN 2158 * 


The particular houndary conditions at y = 0 restrict this development 
Vq 

to the case — =0. 

U 

From the shear distribution (equation (15)) the calculation may be 
made of the coefficients F and S. The atten^t to csQ-culate F and S 
by using the standard profiles together with equations (9) and (I 5 ) was, 
however, unsuccessful for two reasons. First, the ratio 6/0, wMch must 
be known, could not be accurately determined from the standard profiles. 
Second, for reasonable values of 6 / 0 , the calculated values of F were 
positive for values of H for which F should be negative. 

The calculation of the part of F Independent of the shear profile 
was then made both for the standard profiles and the power profiles by 
making use of the kinetic-energy equations (equations (9) and ( 13 )^ 
respectively)} the comparison is shown in figure 3* The closeness of the 
results su gg ests that it is permissible to tise power profiles as an 
approximation for calculating F and S. From equations ( 13 ) and (l5)y 


and 


= -H(3H - 1)[h - 1 - ~" 77 ^ ( 16 ) 

L (H + 5 )(H + T)(H + 9lJ 


S = (3H - 1)[h 2to(3H - 1) ~l 

L (H + 5)(H + 7)(H + 9)J 


(IT) 


The functions F and S, given by equations ( 16 ) and (17)^ respectively, 
are shown in figure 4. 

The fact that the equation 


0 ^ = Rd + S0 
dx 

where F and S are obtained from equations ( 16 ) and (17) ^ respectively, 
does not predict the behavior of the turbulent boundary layer is shown as 

follows: Let cd = 0} then, for H greater than approximately 1.5^ — 

dx 

should be negative. Because S from equation ( 17 ) is positive for 

* 

H > 1.2, it follows that ' is positive. This conclusion is incorrect} 
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therefore, the function for S (equation (l?)) is inconsistent with the 
known behavior of turbulent boundary layers. 

To show that the function for P (equation (16)) is inconsistent 
with the known behavior of turbulent boundary layers, let H *** 1 . 4 . By 

makin g ^ positive and large, — becomes positive and large because 
dx dx 

P given by equation (16) is positive. For positive valTies of how- 

ever, it is known that should be negative. The function for P 

(equation (16) ) is therefore inconsistent with the known behavior of 
turbulent boundary layers. 

In order to detexmiine whether functions for P and S that do not 
result in obviously incorrect conclusions can be obtained“by making the 
shear polynomial satisfy a greater number of boundary conditions at the 
outer edge of the boundary layer, the shear polynomial is generalized by 
writing 


g = (1 - O^ll + a5 + + 1^(1 . ^)A(1 + AO 

2 J 


( 18 ) 


The boundary conditions at the sxirface that are satisfied by equation (I8) 
are 


g = 1 

or T = Tq 

- X 

or = — = 

oy dx 

^ = 0 

or — = 0 


ay2 


At y = 6, the conditions that eire satisfied are 

g = 0 
0 

- 
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In order to evaluate the integral 

^cr 

term — is written as 


r ^ 

Jo 


in equation ( 13 )^ "the 


I = El - X82 


vhere 


81 


= (1 - + A(A + 1){| - A(1 - H 


vA-1 


and, for A ^ 1, . 

gg = (1 - 5)A(l + 2Ai) - 5a(1 + A0(1 - iY 

By using the expression for g^, the equation obtained for S is 

5 f 2 p + 1 Alp r 

(3p + 1 (A + p)(A + p - 1)(A + p - 2) . . . (p + 1) L 


S = 2(3P + i)n 


A(A + l)(p + 1) 


( 19 ) 


A + 1 + p 2(A + 2 + p)(A + 1 + p)J_ 

By using the expression for g 2 the equation for P is found to be 

P = -2(2p + 1)(3p + 1)^ 


K . 

\3p + 1 


A!(p + 1) 


(A + p)(A + p - 1)(A + p - 2) . 


r_i_ 

. (p + i)[a + 1 + 


A(p + 1) 


(A + 2 + p)(A + 1 + p)_ 


( 20 ) 


Hi Pi 
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To avoid positive values for S obtained ffom equation (19) for H < 3j 
it is found that A must he 1 in the expression for It is also 

found that to avoid positive values for P in equation (20) for H > 0, 
A must he «> in the earpression for The values fo;r S and P 

then hecome 


o _ (3H - 1)(H - 1)(H - 3) 

® hTs 


( 21 ) 


P = -H(H - 1)(3H - 1) 


( 22 ) 


The expression for P (equation (22)) is the same as the coefficient of 

— — in equation (13)» letting A — > » makes the coefficient of X in 
U dx 

equation (l8) hecome zero. The shear profile then contributes nothing 

to the coefficient of — — in equation ( 13 ). 

U dx 

Equations (21) and (22) for S and P, respectively, were tested 
by making a computation of H and 0 for the pressure .distribution 

given in table I of reference 10. The computation began at — = 0.075 

c 

with the values given in table I of reference 10. The equations used are 

0 S = -h(H - 1)(3H - l)cD + 0 

dx ' / H + 5 

and 

^ = -(H + 2)cd + 

dx • r 

The equation for — is the Von Karm^ momentum equation. The equation 

dx 

for 0 was obtained from reference 17 and is 

d = o».oo6^3^ 

The cadculated distribution of. H along x was far from the experimental 
curve. 

In an attempt to reduce the sensitivity of the equation for 0 — 

dx 

to the shear distribution, the moment of momentum equation (equation (1^)) 
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in which the shear appears in the coefficient of 0 only as a mean 
value, is used. When the generalized expression (equation (I8)) is used 
for the shear distribution g the result obtained is 


dx 


r 


H(H - 1)(H + 1)^ ^ 3H(H + 1)3 
I 2 (A + 2)(A + 3) 


0 ) + 


(h2 - 1) 






(23) 


where ilr = 0 as required by equation (I8). To keep the coefficient of oi 
negative for all positive values of E, A must equal «> in the coeffi- 
cient of (0. The shear distribution is then independent of the pressure 
gradient. To make the coefficient of 0 negative for values of H 
near 3# A must have the smallest value that it can. take; therefore, 
let A = 1 in the coefficient of 0. ' Equation (23) then becomes 


m ^ -H(e 2 - 1)(H + 1)^ ^ (H^ - 1)(H - 3) gf 
dx 2 5 ^ 


(2k) 


A calcTolation for the example in table I of reference 10 with equa- 
tion (24) resulted in a computed curve for H that was far from the 
experimental curve. 


(b) AssTmiption of constant shear across the boundary layer .- All the 
computations of H have led to values of H much larger than the experi- 
mental values. Therefore, in order to reduce the calculated values of H 
it is necessary to increase S. In order to increase S, the assumption 
of constant shear across the boundary layer is made. For constant shear 
it can be shown that 



by letting g = (l - 0® and taking the limit of the integral as 
B— > 0. Equation (l4), after the assunption of constant shear is intro- 
duced, becomes 


e 


— 

dx 


. (h2 . 1)0 + 


(h 2 - l)i|f 


2 
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In order to make — =0 at H = 1.286 for — = 0, the coefficient 

dx dx ' 

of 0 was arbitr^ily changed to - 1.2862, The equation then 
becomes ^ 

0 — = -H(H + 1)(h 2 - 1) ^ _ (g2 _ 1 . 653)0 + (h 2 - l)ilr ' (25) 
dx 2 

This equation was used for the computation of H with i|r = 0, and the 
results for the example given in table I of reference 10 are shown in 
figure 5 • 

The assomiption of constant shear across the boundary layer was also 
combined with the kinetic-energy equation. When the power profiles and 
the assunption of constant shear are used in equation ( 13 ) ^ the kinetic- 
energy equation becomes 

e ^ = -H(H - 1)(3H - Dm - (H -.1)(3H.-.1J|Z( + (H_+ l)(3H - 1)^ 


The function -H(H - 1)(3H - 1) is shown in figure 3 and the function 

l i ? ■ 1. j in figure 6. When the standard profiles are substituted 

2 

for the power profiles and the assumption of constant shear is made, the 
kinetic-energy equation (equation (9)) becomes 



where the function 

K' 

K • 2 

^ — is shown in figure 6. 


K' 


is shown in figure 3 and the function 
The results of these calculations of 


H 


(with ijf = O) are shown in figure 7* In this case, the use of power 
profiles makes the resvilt somewhat different from that obtained by using 
the standard profiles. 
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would Tae equal to the coefficients of a, 0, and ilr in equation (l4) . 
The ratio of the coefficients of oi is 


H(H + l)(Hg - 1°) 

2 _ (H + 1)^ 


-H(H - 1)(3H - 1) 2(3H - 1) 


The curve of ^ ^ . is given in figure 8 and is seen to be close to 

2(3H - 1) 

xinity. 

The ratio of the coefficients of ■'I' is 

- 1) ^ 4(H - 1) 

(H + 1)(3H - 1) 3H - 1 

- 1) 

The curve of — - — - is also given in figufe 8. The values are far 

3H - 1 

from unity for small valties of H but became equal to unity for H = 3* 
Equating the coefficients of 0 results in 


(h 2 - i)|i + (H + 1) || d^ 


= (3H - 1) 


5 + 




or 


1 1 — 

(h 2 - 1)(H + 1) r 5 || d^ - r ^ 2 |g ^ g2(3 _ 2 ) ( 26 ) 

Jq os 2 Jq 


Now let 


I 


5 ^ dS = a + bH + ch2 

0 
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and 



H -1 

^ ^ ^ dC = J + (^ + ZH2 


T'Hien the shear distribution Is assumed to depend only on the Integrals 
in equation (26) are functions of H alone. Because equation (26) is 
then an identity, the coefficients of the various powers of H can he 
equated to zero. The resulting equations are: 


For 


-a 


1 

2 


= 0 


for 


for h2 


-a-h + 3J-|=0 
a.-b.c-|j + 34-| = 3 


for h 3 
for 

for h 5 


a + h 


c - I q + 31 = -1 


h + c 



c = 0 


The results obtained are: 


a 


28_ 

128 


180 

128 


c = 0 
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i = -22- 

W -I 


51 

128 


Therefore, 


q = _ 


128 


I = - 


ii2_ 

128 


I 


t ^ dt = T - >m 
0 a? 32 


and 


X 


1 lli 


Ilde=-, 


= 7 + ^ 

l6 


Equation (lU), the moment of momentum equation, hecomes 

e = H(bQ - 1)(H + 1)^ . (H - 1)(3H - 1)(T + 22H + 15H2)^ ^ 
djc 2 32 ^ 


(H^ - l)ilr 


(27) 


and equation (13), the energy equation, hecames 

e g = -H(H - 1)(3H - Do. - (H - 1)(3H - 1)(7 ^ 22H ^ ^ 

(H + l)(3H - l)„. 


(28) 


The variation of H with x for the Initial values and the pressure 
distribution given In table I of reference 10 was conqputed by tislng a 

modified form of equation (28).- In order -Uiat — = 0 at values of H 

dx 

In agreement with experiment when o> = 0, the coefficient of 0 In 
equation (28) was replaced by 

(H - Hq)( 3H - Ho)(7 + 22H + 15H2) 


32 
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where 


~ Ho(H0) 

The variation of Hq with Eg was calciilated from the equation 

( 

loSlO ^o = 0*5990 - 0.1980 logio Re - 0 . 0189 (logio 

which was derived to represent a fsdred curve through the experimental 
data (see fig. 9 ) 5 the data were obtained from reference 13 and ftcm 
British results titiat are not generally available. The restilt of a confu- 
tation of H for i|r = 0 and with equation (28) modified as follows 

e S = -H(H - 1)(3H - l)m - - Hq)( 3H - Hp)(7 + 22E + + 


(H H- 1)(3E - 1)^ (gg) 

Il- 
ls given in figure 10, 

Assumptions (b) and (c) lead to somewhat better results than assump- 
tion (a) although they are still not as satisfactory as those obtained 
from the pxirely enfirical relations introduced in references 10 and 12. 

It is clear that this difference is caused partly by the inaccuracy of the 
sinf le assumptions about the shear distribution and can be inqoroved by 
using better descriptions. However, in view of the limited present know- 
ledge of the shear distribution, it does not seem worth while to make 
more complicated assunftlons. 

It may be noted that the final equations obtained for the change of 
the form parameter by the three assumptions are all of the form 


0 g = p(H)(o + S(H)^ 


where a> = — — and 0 = -^. This form is used in reference 12, but a 

TT ^ 


U dx 


pU^ 


different form is used in reference 10, 
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Investigation of Energy Equation 
Since none of the three asstaaptions for the shear distribution ^ 

aTT 

results in. a dependable equation for an investigation is made to 

determine whether a result common to the three assimptlons, namely that 

the coefficient of Tq/pU^ in the equations for S is a function of 

H alone, is very far from true by using experimental data and the 
kLnetic-energy equation without any assumption for the shear. 

If no assumptions other than the boundary-layer assumptions are 
made anfl if in equation (6) n 0 and m = 1, the result is 



If the assumption of a sin^e-parameter family of curves is made 
(f = f(q,H)), then K = K(H), and equation (30) becomes 


e 


dx 



+ CD K + 

K’ 



or, for ijr = 0, 


e g = iZik(H)(l + lo) 


(31^ 


where 

k(H) 

i 

lo 


-(H - 1)K 
K' 


0) 

= T 
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If the assumption Is made that 
f dg is a function of H 
then hecomes 



g = g(n,H), then 
only. Therefore, 


f = f(g,H) and 

5o = ^o(H). Equation (3I) 


0 g = JZfk(H) g + lo(H)] (32) 

• 2 

In order to obtain an estimate of the quantity 1 / f dg \mder 

K Jq 

the assumption that f = f(g,H), reference 12 is used. Equation (7) of 
reference 12 may be ■written as 


0 ^ = JZlki(E)(j - 2.065(1 - 1.4)] 


(33) 


where 


ki(H) = e5(S-1.4) 

Note that Gamer's equation (equation (33)) has the form the kinetic- 
energy equation takes when the assumptions that f = f(Ti,H) and that 
g = g(Ti,l) are used in the kinetic-energy equation. The kinetic-energy 
equation (equation (31)) can also be placed in the form of equation (32) 

when the more general assun5>tion that S = ^ ^'o('n^2) + Fi(n,l) is made 

for the shear distribution. For the purpose of obtaining an estimate of 

the value of 1 - — / f dg, tbe quanti-ty | + In® ^ equation (32) 

K Jo 

is ass'umed to be identical -with the q-uantity | - 2.065(H - 1.4) in 
equation (33) • Then 

|q(H) = -2.065(1 - 1.4) 

and for 1 = 1.5, for example, 

6q(H) = -0.2065 

therefore. 
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or 



f dg = -0.1032 


Therefore 1 



is the difference between two quantities ^ each 


of which is much larger than their difference. It follows that, in order 
to determine values of H not close to separation with any 

accuracy, f and g must he Imown with relatively good accuracy. 


It may he noted that the moment of momentum equation is also sensitive 
to g. This sensitivity can he seen hy writing the coefficient of 0 in 
equation (l 4 ) as 


(h2 - 1) 

/ 


(H + l)(l - 

.7 gat) -1 


\ 



When it is noted that the integral 


■/ 

^0 


g d^- is of the order 


and that H lies between 1.2 and 2.6, the sensltivUy of the 
of 0 to g becomes clear. 


of unity 
coefficient 


In an attempt to determine whether Iq is determined mainly hy H, 
all the data that were used in reference 10 were used to compute Iq hy 
making use of equation (32) in the form 




0 ^ 
dx _ I 
0k(H) ' 


The surface-friction coefficient 0 was calculated hy the formula 
(from reference 1?) 

(k ^ 0.006^33 

and k(H) was calctilated hy the expression obtained from the moment of 
momentum equation 

k(H) = - 1)(H t 1) 


The values of & plotted against H are given in figure 11 . The effort 
to determine whether Iq function mainly of H is Inconclusive. 
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arr 

At least part of the scatter occurs "because ^ 

dx 

from curves -faired throu^ experimental points. 

culation of Iq requires the subtraction of I 
tion which further decreases the accuracy of the 


DISCUSSION 


Althovi^ equation (6) is valid whenever the boundary-layer assump- 
tions are valid, the equations for ^ that result after additional 

assunq)tlons are made do not lead to good agreement with experiment. The 
first of the additional assumptions made is that all velocity profiles 
of the turbulent boundary layer belong to a single -parameter family of 
curves. The experimental data of references J, 10, 11, l4, and 15 
substantiate this assunptlon. 

The second assumption is that the single-parameter family of curves 
can be approximated by power profiles. The data in figure 1, in which 
velocity profiles are compared, and eQ.so the data in figures 2 and 3^ in 
which K and P are conpared, show this assumption to be good, at least 
for H < 1.8. 

From the data in figures 1 to i"fc is inferred that power profiles 
can be substituted for the standard velocity profiles without greatly 
affecting the calculated distribution of H against x for H <1.8. 

To test this inference, the kinetic-energy equation was used with the 
assumption of constant shear across the boundary layerj the result is 
shown in figure 7« As expected from the data of figures 1 to 3^ the 
effect of the substitution of power profiles for the standard profiles 
is noticeable only for H > 1.8. It thus appears that the inaccuracy of 

ATT 

the equations for — that were tested is caused mainly by the surface- 

dx 

friction law that was used and by the assumed shear distributions rather 
than by the use of the power profiles. 

The data of 'references 12 and 15 show skin frictions that increase 
strongly in the region upstream of the separation point before dropping 
to zero- at the separation point. On the other hand, the skin-friction 
data presented in reference lU indicate that the skin friction falls 
monotonically to zero as the separation point is reached. In the present 
analysis a skin-friction law obtained from experiments on flat plates is 
used. It is therefore probable that part of the inaccuracy in the 


and were obtained 
dx 

In addition, the caJ.- 
9 ^ 

from an opera- 

M^) 

calculated values of Sq. 
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equations used to •calculate H is caused ty the use of a relation for 
the skin friction that does not give correct values when there are 
pressure gradients along the surface. 

The assumptions for the shear distribution that were made to obtain 

ary 

a specific equation for ^ were 


(a) The shear distribution depends only on the ratio of the pressxure 


8 

gradient to the skin friction or 

Tq dx 




(b) The shear is constant across the boundary layer 

(c) The shear distribution depends only on the form parameter of 

the velocity distribution ( ' 

Because none of these simple assumptions is derived from a knowledge of 
the details of the turbulent flow, it is not likely that any of them are 
valid. When it is recalled that the coefficient of 0 in both the 
kinetic-energy and the moment of momentum equations is sensitive to the 
shear distribution, it is not svirprlsing that a reliable equation for 

— was not foxmd. 
dx 


In order to obtain a reliable equation for “ flrom eqjuatlon ( 8 ) 

it thus seems necessary to calculate the surface shear and the shear 
distribution across the boundary layer more accurately than in the 
present analysis. Efforts should therefore be made to understand the 
mechanics of turbxilent shear flow sufficiently well to provide an inde- 
pendent relation for the shearing stress that will predict the behavior 
of turbulent boundary layers when used with the Erandtl boundary-layer 
equation (equation (l)). 


■ CONCEUDHSG EEMAEKS 


A general integral form of the boundary-layer equation is derived 
from the Prandtl partial-differential boundary-layer equation. The 
general -integral equation, valid for either laminar or turbulent incom- 
pressible boundary-layer flow, contains the Von Karmm momentum equation, 
the kinetic-energy equation, and the Loltsianskll equation as special 
cases. 
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In an attempt to obtain a practical method for the calctilation of 
the development of the turbulent boundary layer, use is made of the 
experimental finding that all the yelocily profiles of the turbulent 
boundary layer form essentially a single -parameter family. The general 
equation is thereby changed to a sirq)ler one from which an equation for 
the space rate of change of the shape parameter of the turbulent boundary 
layer can be obtained. 

The resulting equation for the space rate of change of the velocity- 
profile parameter is restricted by the assumption that the velocity 
profiles of the turbulent boundary layer can be approximated by power 
profiles. Two of the resulting equations are used to calculate the dis- 
tribution of the profile shape parameter over an airfoil for one experi- 
mentally determined pressure distribution. Although different assumptions 
were tried for the shearing stress across the boundary layer, the cal- 
culated distribution of the profile shape parameter did not agree exactly 
with the e:q)erimental distribution. 

An examination is made of the effect of using the experimentally 
determined single-parameter family of velocity profiles Instead of the 
power profiles on certain functions that occur in the equation for the 
spare rate of change of the velocity-profile parameter. One calculation 
of the distribution of the profile shape parameter over an airfoil is 
also made for the experimentally deteimlned pressure distribution by 
using the single -parameter family of velocity profiles found from experi- 
ment. A comparison of the results with those of a calculation made with 
the same assumptions except for the use of power profiles shows some 
difference near the separation point. It is believed, however, that the 
apparent lack of reliability of the specific equations used to make the 
calculations is catised mainly by the lack of precise knowledge concerning 
the surface shear and the distribution of the shearing stress across the 
turbulent bomdary layer. The present analysis emphasizes the need for 
information concerning the shearing stresses in turbulent boundary layers. 


Langley AeronaTiticeQ. Laboratory 

National Advisoiy Committee for Aeronautics 

Langley Air Force Base, Va., May 22, 1950 
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APPENDIX A 


DETAIIED development of EQUATION (j-) 

Equation (2) can be -written so that terms of the form 
appear e3q>licitlyj therefore, each term -will -vanish at the outer edge of 
the boundary layer. The resulting equation is 

. d“+i)u + + -i- 5:^ + 

m+l|jx' ' 8y J m + 1 dx 

1 ^ (^1 _ ^ ^ uU^+^ ^ ^ 

m + lSy m + 1 Tq dx m + lrodx 

tPu ^ + ^ ^ (Al) 

dx p ^ 


or, af-ter simplification, 

. _L_ri(D*^ - - 

m + iLSx' Sy 

» 

1 (-pntl , dTp _ / mjj _ ^ ^ (J^2) 

m + 1 Tq dx dx p ^ 


Equation (A2) is now multiplied throTigh by y“ and integrated with 
respect to y from y = 0 to y = 6. The resulting equation is 


.^f^yn . u*^)u ay - 

m + 1 uQ Sx 

^ i ^2 r yn(u”i+l - x^^+l)u 

m + 1 ro dx Jq 


_1_ ay - 

m + 1 ^0 ^ 

dy - ^ y“(u™U - utf^)dy = 


1 

P 



ay 

dy 
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or, after simplification and substitution of the formula for the 
differentiation of a definite integral 

^ ay 

the following equation results: 

- .-Tt ■ sris • - .-tijC’" 

.-H s* r • (r*]s ^ i jf fer -!>"«■ 


- / 
pJq 


tFyn dy 

dy 


(A3) 


By integration hy parts, 
^8 


^ yn . t^a+l)v ay = -n jT dy 


and equation (A3) hecomes 
^8 


- m-ti I - iT"]^ ^ X ^ - 

T 4 - (rli ^ s/Br - ti^ ^ = . 


PJo ^ 


dy 


(Al^) 


The velocity v can he eliminated from the term 


r (um+1 _ ^+1^^-1 ^ 

^ r\ 


hy the following development; 
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The velocity v may he written as 

dy + VQ 


_ av 
"" Jo ^ 


or, hy use of the equation of continuity, 

v = -/ ;^dy — / udy + VQ 

Jo 8x To ^ Jo 


, ay - y SJ ^ i 52 . U)dy - « ^ y * VO 

Jo 8x dx To dx Jo To dx 

= /^[u + S(1 - “^ay * y ffi + 2. 52 r/l . HW . 

Jo [_ 8x dx\ uy dx ro dx Jo \ W 


£. ^ 

To dx 


y + VO 


oy a(i - g) ^ ^ ry _ . y(ffl + i 

Jo 8x dxJo I r^dx/ 


or, after terms are collected and. f is substituted, for u/U, the 
result is 

VVU r ay . (g ^ ^ -f r (1 - fiay - yl . VQ 

'-'0 8x \dx ro ^/Uo J 
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The term f vy^"^ dy can now he written as 

Jq 

(ip+l _ unH-l)vyn-l dy = /^ U dy 

Jq ^ ' '-'0 *^0 8x 


'dU U_ ^o \ 

dx ro dx 


) Jq ~ ^ 


+ Vq 


^dy 


or 


. um+l)^-l ay = jf^2 yn-1^1 . fm+1) T) d^dy + 

(i ^ I [f<" - ^ 


Torf“-*-^y'V'^(i - f"V 


Now let 




8(1 - f) 
Sx 


dy 


dy = Ii0 


n+1 


r y“-^(l - /^(l - f)'3y 

Jo Jq 


dy = J0^+^ 


J j^(l - dy = M0“+^ 
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The term dy now become s 

& 

J' (um+l _ ^m+l)vyn-l dy = + IL + 




Now let 


r y“(l - f*a+l)f dy = 
^0 


m 


n+l 


and. 


08 

/ (f - f“)yn dy = L0“+^ 
Equation (a 4) can then be written as 

_ 1 + s 

m + 1 dx' ■' m + 1 ■*■ 


m 


m 


jSL-(^ + 2L - 11)0’^+^ + — S~ Vntf^^Q0“ - 

+ lldx Tq dx / m + 1 y- 

1_ ^2 + iT-^ ffi L9“*^ = r® ffyn L’-ay (A5) 

+ 1 To dx dx P '-'0 ^ 
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in 


After — Is escpanded and terms are collected, equation (A5) 

dx 

becomes 

(n + 1)N — + ef— - nil) + - — lN(m + 2) - n(j - M) - 
dx \dx -V U dx*- 

L(m +1)] + — ^|w - n(J - M)] - nQ ^ = 

-* rQ dx •- u 

^ p6/e 

-(m + 1 ) -0- / f“ijn ^ dTi 

pU^ ^0 

where 


on 


n =3L 

e 
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APPENDIX B 


SIMPLIFICATION OF TERMS IN EQUATION (8) FOR POWER PROFILES 


Calculation of — + nl 

an 


. The definition of is 

Ne °-+1 = dy 

therefore, 

^ 8/0 


N = y' (1 - f^l)ftiJi dTi 


and 




or 




The definition of I©^'*'^ is 


10 


n+1 




therefore. 





KACA TU 2158 


^3 


Then 


~ + nI=r^ ^E-(m + 2)f”^+^ nil dTi + 
SR 


(bi) 


By integration hy parts, 
5/0 


n I an = -jyii . - 


(m + l)f“ 


^ f ^ dnl 

^ '-'0 J 


r\°- dij 


(B2) 


When equation (B2) is substituted into equation (Bl), the following 
equation is obtained; 

5 C' lo I " i)^" 


Use is now made of the power-law assumption 


f =5^> = nP(| 

^O, 


0\P 


Then 


1 - 2p^ ^ P ^ ^Plog i 


Bf 


an 2(p + l)(2p + 1) 


and 


^ ^ 
^ Jo 


an - f g 


1 - 2p^= 


1 - 2p^ ^2p ^ g'^log i 

2(p + l)(2p + 1) 2 _ 


2(p + l)2(2p + 1) 2[p + 1 (p + 1)^ 




V - 



After a lengthy manipulation, ^ + nl is found to be 


— + nl = ^ (p + l)°^(2p + l)^(pm + n) 


as 


p’^'*’^ jp(m + 2) + n + 
where use has 'been made of the following equation: 


S = (p + l)(2p -»• 1) 
e p 


Calculation of N 
The definition of is 

.8 

'0 


= r (l - fy^ dy 

Jo 


When f = ^ is used, the equation for N is 

jj ^ P(m + 1) 

jp(m + 2) + n + l| (p + n + 1) 


or 


N = •»• 1)(P + 1)°^^ 

p!*^[p(m + 2) + n + ^{p + n + 1) 


Calculation of J 
From the definition of 


J0 


n+1 


= (l - y\l - f)d^y^-l dy 
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\Then f = the equation for J is 




.or, after a lengthy manipulation. 


_ (p + l)n(2p + 1)11+1^^^^^^ (p+ 1 +n) |^+n+ 2 (p+l)] +p(m+l)+l + n 

(n + l)(p + 1+ n) jp(m + l) +"n + ^ |p(m + l) + p+l + ^ 


Calculation of M 


From the definition of 

pB 

= / (1 - f“+l) yn dy 
% 

VThen f = ^P, the equation for M is 

. M= (m+l) ■ (P 1)*^^ 

p“(n + l)|p(m + 1 ) + n + :^ 


Calculation of L 


From the definition of 


- f”‘-l)fyii dy 


\Jhen f = ^P, the equation for L is 


L = (m - 1) (P l)°-^^(ap -■ 

p^(p + n + l)(pm + n + 1 ) 



Calculation of 


H(H - m) - n(J - M)(H + l) + (m + 1)L 
dK 


dE 


+ nl 


Frcm tbfi e:xpre salons for J and. M, the expression for J - M Is 

j _ M = (P + 1)^(2P + (p + 1 + n)[pn + n + 2(p + lH + Gp(m + l) + 1 + n] 

(n + l)(p + 1 + n) |p(m + l) + n + ]J [p(m + 1) + p + 1 + n] 

* 

f ^ , X (p + l)^-*-^(2p + l)^-*-^ 

pJ^(n + 1) [p(m + 1) + n + ^ 

After a lengtlny simplification, the result Is 

-(m + l)(p + l)'^(2p + 


J - M 


pP(p + 1 + n) ^(m + 1) + p + l+ i^ 


or 


-n( J - M) (H + 1) 


2n(in + l)(p + l)^'*’^(2p + l)’^'*’^ 


p*i(p + 1 + n) [p(m + 1) + p + l + ry 
where H 2p + 1 was used. The esqsresslon obtained for N(E - m) + (m + l)L Is 

(m. + l)(p + l)°*^(2p + l)^'*'‘*^(pni + m + n) 


N(H - m) + (m + 1)L = 


- \ W— M 


^ip + n + i; |p(m + 2) + n + (pm + n + i; 
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and the eacpreaaion obtained for ir(H - m) - n(J - M)(H + l) + (m + l)L la 

»(H - m) - n(J - M)(S + 1) + (m + l)L = &(l- n -H) H- p(p. I- m t n)] 

P^(p + 1 + n) [p(m + 2) + n + ^ (pm + n + 1) 

By Buhstltution and simplification 

g(H - m) - n(J - M)(H + l) + (m + 1)L ^ -4p(p + l)(2p + l)fpCm +2) + n + j 

. pm + n + 1 

dH 


It can alao be shown that 


r®/® S 

+ n(j - M) - (m + 1 ) °g di] ^p 


dH 

dH + 


pm 


^|p(m + 2) + n + i] (ap + 1 + 


[p(m + a) + n + Ci“«| 5 a£ 


Evaluation of -N + n(J" - M) + nQ 
Prom the results for N and J - M 


-N + n(J - M) = ±. 

p*^jp(m +2) + n + ^ 


*r 
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For Q, the developnent is: 


/-«5 

qp^= (1 - 

dn 


cly 


anfl with f = the following expression is obtained for n ^ 0: 

Q = 1 (P + 1)^(2P + 1)^ 

“ p’^"^[p(m + 1) + 

Then, by siibstitution and simplification, for n ^ 0, 

-(m + l)(p + l)“^^^(2p + l)’^(ian + n) 


-N + n(J - M) + nQ ‘= 


p'^jp^m + 2) + n + Jjp(m + 1) + 


(B3) 


If use is made of the previously derived result that nQ = 1 for n = 0, 
the following equation is obtained for n = 0: 


-K + n(J - M) + nQ = t 

p(m + 2) + 1 


(B4) 


If n is placed equal to zero in equation (B3), equation (b 4) results; 
therefore, equation (B3) is valid for n = 0 as well as n 0. 

Then, for all velties of n, 

-N + n(J - M) + nQ ^ 2(p + l)p[p(m + 2) + n + j 
dU 


p(m + 1 ) + n 
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Figure b.- Values of S from klnetlc-euergy equation vlth the assumption 

of constant shear. For power profiles, S = - */.• for 

K •• 2 ^ 

experimental profiles, S « - — ~ — . 






Figure 7 .- Comparison of experimental values of H and values of H 
calculated 1)7 Idnetic-energy equation with constant shear across 
boundary layer. 
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Figure 8,-. Satlos of pressure and suction factors in energy and moment of 

momentum equations. 
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Figure 9«“ Variation of with Eg. 













